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I. INTRODUCTION 



1. The basis of relativistic superfluid mechanics was emphasized by Israel |T[| and Dixon 
and their achievements were developed by several authors p-^H,?^. Since the rela- 
tivistic superfluid matter of neutron stars, being the main object of applications of this 
theory, contains quantum vortices, an ample discussion of a relativistic superfluid vortex 
prompted in diverse forms P-jUfl and taken into account by Carter and Langlois [H] in 
a new method designed for general hydrodynamic description of a relativistic superfluid 
with quantum vortices P~5|JTS|JT1^JTT[ | . For the previous discussion did not explicitly include 
vorticity whose presence changes the equations of motion [14|. The latter model shall be 
applied in the present study to investigation of the sound propagation through a superfluid 
containing quantum vortices. It should be noted that the sound |L9],[7|,[13| and shock wave 
propagation |20| in a two-constituent relativistic superfluid have been already discussed (as 
well as the recent analysis of the Cauchy problem [?]); the sound in a "vorticized" relativistic 
superfluid, however close it seems to the usual first-second sound, bears principal difference 
from the former, because vorticity w vp is interpreted as a tensor dynamical variable on which 
the Lagrangian function depends | 15f| . So, in the model which we shall exploit below |T5|JT8[| 
the role of the second constituent is played by the vorticity considered as a vortex " liquid" 
rather than a vortex lattice where the so-called Tkachenko-Dyson waves |2T| can be observed 
(they are, evidently, inconceivable in the present model). Thus, we consider the waves in a 
relativistic superfluid in the "continuous" or hydrodynamic limit. 

We use a natural system of units Ti = c = 1 and a metric g vp =diag( — 1, 1, 1, 1). 



II. RELATIVISTIC SUPERFLUID DYNAMICS WITH QUANTUM VORTICES 



The equations of a relativistic superfluid with quantum vortices |L8[ include the equation 
of motion 



f w va = 



(1) 



2 



and the conservation law 

V u n u = (2) 
of the particle number current n u . The total current 

f = n » + V a \ va (3) 



is also conserved: 



v„r = o 



The vorticity 2-form 



w vp = 2V[„ /If,] (4) 

satisfies the closure condition (square brackets imply antisymmetrization of indices) 

V [Q w up] = (5) 

but it is not trivially zero, for the vortices are presented. 

Since the current n u and the vorticity w up are considered as new dynamic variables, a 



variation of the Lagrangian L is given by formula flTJ 



dL = fi u dn" — -X up dw up (6) 
Instead of the proper Lagrangian L we can use the pressure function 

* = L - ii v n v (7) 

whose variation 

d^ = -n v dii v - \\ vp dw up (8) 

is expressed in terms of new primary variables \i v and w up . In general the pressure depends 
on three invariants 



I'- 



-w up w 



up 



h 2 = h v h, 



(9) 



with helicity vector h u , whose conservation 



V„ h v = 



(10) 



is a consequence of Eq. (f|) and fl5|). The secondary variables n u and are expressed 
through primary variables by relation 



\ ( „ \ 



n 



with 



<9^ 

pva _ 2 gUa 



<9* 

oh 2 



R ai/p G U ^ P1 



/'a 



(11) 



(12) 



(13) 



It should be emphasized that the vorticity 2-form does not reflect the fine structure of 
the vortex cell, since the system of vortices is considered as a vortex " liquid" . All variables 
are averaged over the vortex cell and vary over the macroscopic range. In the non-relativistic 
limit only space components of the vorticity form w up survive, while uvo = 0. 

III. INFINITESIMAL DISCONTINUITIES 



Let be a space-like normal to the front of infinitesimal discontinuity. We shall apply 
the so-called Hadamard technique |22| for investigation propagation of small-amplitude per- 
turbations. Recently Carter and Langlois |7j succeeded in calculating the first and second 
sound speed in a relativistic two-constituent superfluid. We use the same method to find 
the sound speed in a superfluid containing quantum vortices. The infinitesimal deviation A 
of an arbitrary quantity A is proportional to the deviation of its gradient 



[V, A] = AT U 
4 



(14) 



Applying formula fll4|) to the conservation laws ( PH) and (0) we get, respectively 



r„ = o 



(15) 



Y v n v = 



(16) 



The closure condition @ gives equation 



Y[ a w up ] = 



(17) 



while the equation of motion (JTJ) and (ffl) in the linear approximation yield 



n v T[ v fi p ] + Y a \ va w vp = 



(18) 



A relationship between the deviations of the primary and secondary variables is estab- 
lished by the formula 



ftupa Qvpr\$ 



obtained by differentiation of equations (|TT|)-(|13|) where 



d 2 ^ d 2 ^ / a n 



d^ a ~d(ji 
d 2 ^ 



d[i 2 dh 2 



(19) 



(20) 



d[i 2 dW 2 



d/2 2 dh 2 
d 2 V 



dh 2 dW 2 



d(h 2 



M 2 



h K e KV1Td 



(21) 



(22) 



<w 2 



~<9(iy 5 



- 4 



<9 2 * 



dW 2 dh? 



<9 2 ^ 



d(fr 2 ) 2 dh j 



(23) 
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Substituting (|I9|)-(|2"3"1) in (p~5|)- (P~8|) we write ten equations 

e ua ^ T u {faw Pl + fi a w M ) = (24) 

^[ a w V p] = (25) 
T v (F™fa + CT^w^) = (26) 

(R vpa fa + C^Wtf) w va = (27) 

with ten independent unknowns, namely fa, fa, fa, fa, w i, u> 2> w 03 , w 12 , w% 3 , w 2 3- A 
linear system (p4|)-(f2"7j) is consistent if and only if its determinant vanishes, that yielding the 
characteristic equation for speed u of the wave. 

IV. DILATONIC MODEL 

In order to obtain an explicit analytic solution we shall confine ourselves with a particular 
dilatonic model |18[ whose Lagrangian L and, hence, the pressure function \l/ do not depend 
on the helicity h. Although this model is devoid of complete strictness, for weak vorticity is 
provided (GW 2 <C fin) , it corresponds to the conditions inside the neutron stars and, hence, 
describes satisfactory the real superfluid matter. Thus, formulae (120)- (EST) reduces to 



d 2 ^ 

r -V""W (28) 

^ = "w^" (29) 

<9vfr f) 2 \tj 

G vprr& = -2——g uri gP» - 4— *w vp w n6 (30) 

dW 2iJ y d{W 2 ) 2 

while ( p2|) remains the same. 

For a dilatonic model it is convenient to practice with a dilaton <3> 2 instead of the set ([5|). 



In the weak vorticity limit |18[ we have 



* [$ (//, W)] = $ : 



d$ 2 



-V[<&] 



Since 



thereby, 



and 



n 2 = 2$ 4^ + 2KW& <$>~ 2 = ^ 



d& 



(i 



n 



} dV 



+ 2KW 



d$ 2 



d » 2 -w dw \ 



/2V' 1 



where 



V 



dV 



(31) 



(32) 



(33) 



(34) 



(35) 



and quantity K (whose logarithmic dependance on vorticity is neglected), being proportional 
to the circulation quantum round an individual vortex, satisfies the inequalities 

KW < V <$> 2 KW <C V" (36) 

while the sound speed is given by the formula 

2 ndn 



$ 2 V" 

" s ~ Jidn ~ $ 2 V" + /i 2 
we obtain the derivatives 



In the light of ([31])- 



dm dm fdfi 



I 



<9/z 2 d<$> 2 ' \d§ 2 



dm dm 



fd$ 2 \ 
dW 2 d$ 2 \dW 2 ) 



and 



dm 



1<9$ 2 



i 



d (/i 2 ) 2 2 dfi 2 AV" 

dm K 
d^ 2 dW 2 ~ ~4WV" 



1 K 



- - ( $ 2 + 
d^ 2 ) 2 4W^ 3 1 " 



A" 



WV" 



(37) 

(38) 
(39) 

(40) 
(41) 
(42) 
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V. EXPLICIT SOLUTION 



Let the vortices be aligned along the axis z and there is no dependence on z. The only 
non-zero components will be u> i, w 02 , ^12- Also only fi and n° are not equal to zero in the 
local reference frame commoving the fluid (or one may merely consider the fluid at rest). 
We can choose the normal vector T v as 

V v = (u, cos x, 0, sin x) = {u, S, 0, fi) (43) 

where \ is the angle between the axis z and the direction of the wave propagation. The 
coefficients (p8[) -([30|) and (p2|) , then, take the form 

2 uO aO 

F ua — ~$ 2 g ua — 9 v f (44) 

Q»** = Eg v0 w r * (45) 



(46) 



W 2 ) 

where H = fj, K/(W V") and G = $ 2 K/W. Note that in Eq. (|46|) we do not add the negligible 
term KW/V" incorporating Eq. ((4"2[). Substituting formulae (|4^-(|4"6|) in equations (|24] ) -(|27 D 
and taking into account (|9|), namely 

W 2 = w 10 w 10 + w 20 w 20 + w 12 w 12 (47) 

we get 

- jjLoVtWu - fllQw 02 + fa^Woi - feU + W 2 3^fi0 + Wi 2 {l[i = (48) 

- uw 12 + SW20 = (49) 
- uw 13 + Sw 30 + ilw i = (50) 
-uw 23 + ttw 02 = (51) 



S 



T,W 23 + &W 12 = 



(52) 



where 



- ufio F 00 + S/ixF 11 + fi/isF 33 - u (Q 001 w 01 + Q 002 w 02 + Q m2 w 12 ) = (53) 

fi Q Ro + <Xq = (54) 

/t (-2n°£ + R x ) - jj^nPu + af = (55) 

/to R% — 2un°fi 2 + Wrft a 2 — (56) 

/t O + /t 3 w = (57) 

a f = Gaf (58) 

#0 = Zw 02 R 120 - uw ok R 0k0 R2 = W02ZR 010 + w 12 uR 0W 

Ry = Zw 01 R 0W + w 12 (ZR 120 - uR 020 ) (59) 

af = -Bq U af = -B U ag 3 = 

a 12 = (1 - B ) U af = -w 01 Q af = -w 02 Q (60) 
af = B X U af = (1 + £1) 17 af = fiw i 

a} 2 = Si £/ a} 3 = af = -Qw 12 (61) 
a^ 1 = (1 - 77 2 ) 17 af = ~B 2 U af = nw 02 

af = -B 2 U af = ttw 12 af = (62) 



while 



w 01 w 12 w m w 02 w 01 w 01 . . 

^ = —^2— B ^ = —^F2— B ^ = —^F2— ( 63 ) 

U = uw 12 + Sw 02 (64) 
It should be noted that R\ in (|55|) , (|j9|) and Q 012 in ( p3l) survive in the non-relativistic limit. 
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VI. WAVES 



A. Horizontal waves 



We call the wave propagating along the axis x (i.e. at the right angle to the axis of 
vortices z) "horizontal" since for this case Q = 0, X = 1. The matrix of the linear system 



(48-57), then, has the form 



t 










—Wo2 — UW12 

















\ 

/^o 

















1 





u 


























1 





u 
































u 





























1 


-uF 00 


F n 








Q 01 


Q 02 


Q 12 








Rq 











«0 01 


af 





< 2 








-1 + Ri 


— u 








< 


af 













R2 





— u 





< 


af 





&? 

















u 

















J 


variable 


we 


introduced the notation = 


0/(2n 


) and Q 



(65) 



solution follows from a condition of the determinant ( |65"D vanish that implies the possibility 
of longitudinal waves (/ii ^ while fi 2 = A3 = 0, also w 03 = w 13 = w 2 3 = 0) with the 
relevant speed found from the characteristic equation 



( u 2 F 00 -F n )U 2 det 



-B 
l-B 2 



1 - B (1 - u 

B 2 (1 - u) 



J 



o{K) 



(66) 



which, according to Eq. (0), gives the usual sound speed (P7|), for 

,2 



^00 



V" 



We neglected in ( pB]) the insufficient term R\ which reflects the negligible space anisotropy 
of the first sound; in the relevant equation for the first sound in vertical direction (see next 
section) this term is absent at all. The second solution of (|66|) is 



10 



20 



w 



(67) 



In the non-relativistic limit ioqi = W02 = and the only non-zero elements of (|60| - |62l) will be 
a? 2 and a\ 2 . The matrix (|65|) , then, determines the only "first sound" solution, while (|67|) 
does not occur, for it pertains to a relativistic superfluid. 



B. Vertical waves 



Let the wave propagating along the axis z of the vortex cell, when £ = and O = 1, let 
be called "vertical". Eq. fl52|), then, gives wi 2 = 0, while the variables i&oi and wq2 can be 
excluded by means of Eqs. fl50|) and (^1|). Then, we have to calculate the determinant of 
matrix 

/ 



/if) 


/'1 


f'2 




^03 


Wi 3 


W23 


-W12 




WQ1 


-UWu 











-uF 00 








F 33 





-uQ 01 


-ug 02 


-1 








— U 











R 














«g 3 + ua^ 1 


a 2 , 3 + Ma 


Ri 


— u 








Gw 01 


ua? 1 


a 23 + ua 


R2 





— it 





Gw 2 


«2 3 + ua® 1 


ua , 2 



02 





(68) 



where the first row includes the independent variables. The speed of the longitudinal wave, 
characterized by fii = fi 2 = 0, follows from the equation 



/ bl b 2 ^ 



(u 2 F°° - F 33 ) det 



"0 

1 1,2 



Wqi b\ b\ 
w 2 b\ b\ 



(69) 



hi 



-Wq\ — uB U 



-W02 - uB U 



b\ = uB x U b\ = -w 12 + u{l + B 1 )U 



(70) 
(71) 
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b\ = w l2 + u (1 - B 2 ) U b\ 



-uB 2 U 



(72) 



and coincides approximately with the usual sound speed fl37|) . For the determinant in ( p9|) 
as one can check by means of ( |63[ ) and (|?0|)-(|72|), does not vanish. 

As for the transversal waves ((13 = 0) they are determined by equation 



/ 



det 



V 



b 1 



b 1 



b 2 



Ir 







(73) 



-w 2 wqi 

-uQ 01 -uQ 02 



— it Wqi b\ b\ 

— it w 2 

which trivially equals zero (due to the matrix degeneracy in (|73|) , implying impossibility of 
transversal waves. 



VII. CONCLUSION 

In order to investigate wave propagation in a relativistic superfluid with quantum vortices 
we have derived the linear system (|2~iD-(PTD by means of the Hadamard method and 
have solved the appropriate characteristic equation. For a dilatonic model in the weak 
vorticity limit we found two types of waves propagating at the right angle to the axis z 
of vortices. Besides the usual sound (|37|), an additional horizontal (in direction orthogonal 
to z) "second-sound" branch (^7]) exists, while only the usual sound may propagate in the 
direction z. It should be noted that no wave propagation with the speed close to K W/V or 
Q 2 K W/V" were found (the latter term determines the anisotropy of the first sound in the 
horizontal and vertical direction; we do not include this term in Eq. ( |66D because it is small 
and does change the result on qualitative level); for the second constituent was treated as 
a vortex "liquid", no transversal waves occur. On the other hand, due to orientation in z 
direction the condition of wave propagation along z and x (or y) is different. However, only 
the usual sound may propagate through a non-relativistic superfluid, and that follows from 
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Eqs. ( f!8[ - |57D taken in the non-relativistic limit (ww 
non-relativistic system 



0). Because the matrix of the relevant 
\ 



Hi) 


Hi 






111m 






llh o 


^ lo 


7/Joq 


— riwi2 


o 





— UW\2 
















S/in 


o 














-S 





— M 








o 











n 





-S 





— It 





o 











o 


Q 











—u 























n 





E 


-uF 00 







QF 33 











-uQ 012 





































+ \w l2 R lw ) 


— u 






















— flw-i 








—un° 





-uw l2 











VIW12 





n 








u 





















V 

yields for the horizontal waves (Q = 0, S = 1) the characteristic equation 

»»*-+^ (1-^ = 

where 



12 



(74) 



(75) 



R 



[20 



in accordance with (B9| 



WV"" ~"2ri° 2W 
1, and merely 



<C 1 



(76) 



U 2 F 00 + F 33 = q 



(77) 



for the vertical waves (Q = 1, S = 0). Indeed, Eq. flT5j ), which we write here in the exact 
form with the negligible term ([76]), follows from (|66|), while Eq. ( [77D follows from fl6"9|) and 
they both determine the usual sound speed fl3~T|). 

Nevertheless, formula QB"7| ) requires more discussion. What if (|B7|) will be of order of the 
usual sound (^3)? Then we have obviously missed the some information by vanishing wq 2 . 
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We do not know the explicit form of W02 and wu', although u>02 vanishes in the non-relativistic 
limit, the ratio flB7|) . The analysis of the exterior region of a vortex line |H] is not sufficient 
to determine the components of the vorticity 2-form. These components are generated by 
the vortex core and they can be derived from its internal structure; because outside the 
vortex core, i.e. outside the support of the vorticity 2-form w up , the irrotationality condition 
w vp = takes place. So, the complete answer requires explicit derivation of u>02 and wu- 
Of course, it seems reasonable to discuss in future the shock waves. However, equations 



(|27|) have already convey [or: reveal/ the absence of any other modes except (j37|) and 
For instance, the waves through the pure vortex constituent (fx v = 0) or those similar 



to the fourth sound |2(J in superfluid (w up = 0) are impossible. This may serve a hint to 
search possible shock wave solutions. Moreover, the analysis of a more complicated model 
with sufficient vorticity W and the pressure function \1/ depending on the cross term h may 
also be performed by means of Eqs. (|24])-(p7|). 
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